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Abstract 

We propose a novel type of Hamiltonians for quantum field theories. They 
are mathematically well-dehned (and in particular, ultraviolet finite) without any 
ultraviolet cut-off such as smearing out the particles over a nonzero radius; rather, 
the particles are assigned radius zero. We describe explicit examples of such 
Hamiltonians. Their definition, which is best expressed in the particle-position 
representation of the wave function, involves a novel type of boundary condition 
on the wave function, which we call an interior-boundary condition. The relevant 
conhguration space is one of a variable number of particles, and the relevant 
boundary consists of the configurations with two or more particles at the same 
location. The interior-boundary condition relates the value (or derivative) of the 
wave function at a boundary point to the value of the wave function at an interior 
point (here, in a sector of configuration space corresponding to a lesser number of 
particles). 
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1 Introduction 

In most quantum field theories (QFTs), the formulas that one obtains for the Hamil¬ 
tonian (by means of quantization or other heuristics) contain terms for the creation 
and annihilation of particles that are mathematically ill-defined; for example, they can 
be written as integrals over all wave vectors fc G but the integral is divergent at 
|fc| —)■ cx) (“ultraviolet (UV) divergent”). The problem can be avoided by discretizing 
space or treating the electron (and other particles) not as a point but smearing it out in¬ 
stead over a small positive radius (“UV regularization” or “UV cut-off”); however, these 
procedures tend to break the Lorentz invariance, and there is no empirical evidence for 
either discrete space or a positive electron radius. 

We describe a novel way of defining a Hamiltonian for a QFT and report that we have 
proven for some model QFTs, in joint work with Jonas Lampart and Julian Schmidt 
PMEU, that such Hamiltonians are well-defined and self-adjoint. In our approach, 
space is continuous (as opposed to a lattice), and the radius of the electron (or other 
particles) is zero. The key element of the approach is a new type of boundary condition 
that has not been considered before in either physics or mathematics, as far as we 
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know; it was first considered by one of us in [TH] for a different purpose. We call this 
condition an interior-boundary condition (IBC) because it relates the values of 'ip on 
the boundary of configuration space Q to the values in the interior of Q, as we will 
explain presently. The IBC is formulated in the particle-position representation of the 
state vector pj in Hilbert space Here, “particle representation” means that is 
represented as a Fock space (or, if appropriate, a tensor product of several Fock spaces), 
and “position representation” that the contribution from the n-particle sector of Fock 
space is represented (like a wave function in quantum mechanics) as a function of n 
points in 3-dimensional physical space. Specifically, if ^ is a single Fock space then 
the Ip G can be viewed as a function on a configuration space of a variable number 
of particles, such as 


OO CXD 

Q = U 2- = U \ 

n=0 n=0 


( 1 ) 


see Figure [T} where 


A 


n 


,X 


i) G (M^)” : Xi = Xj for some i 7 ^ j| 


( 2 ) 


is the “diagonal,” i.e., the set of collision configurations (i.e., those with two or more 
particles in the same location). The relevant boundary dQ of Q is U^qA„; the IBC 
relates the values of xp on dQn = A„ to the values of tp in the interior of Qn-i, namely 
at the configuration with one particle removed (or possibly with more than one particle 
removed, if more than two particles collide). 

Our approach, which we call the IBC approach, provides a Hamiltonian Hjbc whose 
domain is defined using an IBC. It is often assumed that the Hamiltonian of a QFT is the 
sum of two self-adjoint operators, the free Hamiltonian and the interaction Hamiltonian. 
In contrast, it is a feature of the IBC approach that its Hamiltonian Hjbc cannot be 
split in this way, as we will explain in detail below. 

Work in progress on Hamiltonians with interior-boundary conditions includes |9l 
ESI [ini [m 13 0 [13. Mathematical proofs are provided in EllIDllII]; we give a gentle 
introduction to IBCs and the ideas behind them in [16] . 

We hope to apply the IBC approach to quantum electrodynamics and other serious 
theories in the future; work on implementing IBCs for the Dirac equation is in progress 
HH. For the time being, we report results for simple non-relativistic model QFTs. We 
describe examples of IBCs and how they help define a Hamiltonian, results about the 
rigorous existence and self-adjointness of the Hamiltonians, and how these Hamiltonians 
are related to some known cases in which a UV cut-off can be removed, thus making 
it plausible that the IBC Hamiltonians are physically relevant and not merely mathe¬ 
matical curiosities. Specifically, we consider two models: In Model 1, x-particles can 
emit and absorb y-particles, and both kinds of particles are non-relativistic; we give the 
full definition in Section below. Model 2 is a simplihed version of Model 1 in which 
the a:-particles cannot move but are fixed at certain locations; this arises as a limiting 
case of Model 1 in which the mass rux of the x-particle tends to 00 ; for simplicity, we 


3 



(b) 


(a) 


(c) 



Figure 1: Illustration of a configuration space for a variable number of particles in 1 
dimension: (a) The 0-particle sector contains only one configuration, the empty conhg- 
uration; (b) the 1-particle sector is a copy of physical space; (c) the 2-particle sector; 
(d) the 3-particle sector. In this example, we have not removed the diagonal A„. 


consider only a single x-particle (and call its location the origin). That is, ^/-particles 
can be created and annihilated at the origin of For Model 2 (the one with hxed x- 
particles), it is known that, after a UV cut-off is introduced, the cut-off can be removed 
by means of a suitable limiting procedure (renormalization), and we can show that the 
limiting Hamiltonian Hoc thus obtained coincides with Hjbc up to addition of a hnite 
constant (see Section below). For Model 1, however, such a procedure is not known 
(although it is for a similar model [H]). We are currently investigating whether the 
IBC approach provides, also for Model 1, a well-dehned and self-adjoint Hamiltonian 
without a renormalization procedure ra. 

This paper is organized as follows. In Section we describe two simple models of 
non-relativistic QFT. In Section]^ we describe for each of the two models how to set up 
an IBC and the corresponding Hamiltonian. In Section we develop a general theory 
of IBCs under the simplifying assumption that the boundary has codimension 1; this 
assumption is actually not satished in the application to QFT without UV divergence, 
but Section 1^ adds to the full picture about IBCs. 
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2 Two Models of Non-Relativistic QFT 

In this section, we describe two simple non-relativistic model QFTs. We set up their 
Hilbert spaces 3^^ describe their original Hamiltonians Hong, and explain why they 
are UV divergent. Model 1 is a natural, physically reasonable model, while Model 2 is 
an artihcially simplihed version that will be useful for an easy discussion of IBCs. In 
Section we will describe, for each of the two models, our alternative definition of the 
Hamiltonian involving an interior-boundary condition. 

2.1 Model 1: x-Particles Emit and Absorb ^-Particles 

Model 1 is a QFT adapted from [Ul p. 339], [H], and the Lee model [12]. It involves 
two species of particles, x and y; the x-particles can emit and absorb ?/-particles. For 
simplicity, we take both species to be spinless and assume the non-relativistic dispersion 
relation E = p^/2m for both the x and the y-particles with masses mx,my > 0. The 
Hilbert space is a tensor product of Fock spaces, 

^ 0 (3) 

with 

OO 

= 0 S±L\R\ C)®” , (4) 

where S- is the anti-symmetrization operator, 5+ is the symmetrization operator, and 
are their ranges (i.e., the spaces of (anti-)symmetric functions on 
Here, we take x-particles to be fermions and ^/-particles to be bosons]^ is understood 
as physical space, i.e., ordinary position space. In the following, we simply write 
for C). A vector T G 3^ can be regarded as a function 

V' • Qx X Qy —)■ C (5) 

with Qx = Qy the configuration space of a variable number of particles, 

OO OO 

S. = S,= [JSn=U(*")”' (®) 

n=0 n=0 

where the union is understood as a disjoint union and (M^)° = {0}. We will discuss 
removing the collision conhgurations later. 

We call the function -0 the particle-position representation of the vector T G 3^. 
A generic element of Qx x Qy can be written as {x,y) = (a^i,..., 2/i, • • ■, 2/n)) 

where bold-face symbols denote vectors in 3-space, while x denotes a configuration 
of x-particles and y one of y-particles; we will often write x”^ instead of x to convey that 

^This choice is contrary to the spin-statistics theorem; but that does not matter for our purposes, 
as the latter pre-supposes Lorentz invariance. 
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the configuration consists of m x-particles, and likewise y" instead of y. We call 
the m-particle sector of Qx and x the {m,n)-particle sector of Qx x Qy. 

Likewise, we say {m,n)-particle sector of (or of \k, or of ip) and write (or 

^(m,n) Qj, 0 S+L'^(MP, C)®”, respectively for the projection of 

'h to that subspace, and for the restriction of ip to (M^)™ x (M^)"'. The ip^^’'^l function 
is anti-symmetric in the x variables and symmetric in the y variables. 

The spaces Qx, Qy are equipped with the volume measure 

OO 

t^{S) = n Qn) for 5 c Q, . (7) 

n=0 

The inner product in is then given by 


{iP\P)= fi{dx) n{dy)ip*{x,y)p{x,y) 
J Qj. J Q. 

OO OO 


y 


HU 


m=0 n=0 


I;2(R3m+3n) 


( 8 ) 

( 9 ) 


2.1.1 Original Hamiltonian 

As our example of a (non-relativistic) Hamiltonian Hong, we take 

fc2 ”9, fc2 ”■ 

^ i=l ^ j=l 

+ nEQip{x'^ ,y'^) 

m 

+ gVnTT'^ip{x^,{y^,Xi)) 

i=l 
m n 

+ ^ ~ \ yj) ’ 

V ^ i=l j=l 

using the notation for the 3-dimensional Dirac delta function and 


y"" \ Vj = ivi, • • •, Vj-i. Vj+i, 


( 11 ) 


for the configuration of n — 1 ^/-particles with the j-th particle removed; Eq is the energy 
that must be expended for creating a ^/-particle (the “?/ rest energy”), and G M is a 
coupling constant (i.e., the “charge” of an x-particle). The first term in (10) is the free 
fermion Hamiltonian Hxip, the second and third terms are the free boson Hamiltonian 
Hyip, and the third and fourth line together are the interaction Hamiltonian Hi^tg^ip 
responsible for the creation and annihilation of ^/-particles. In terms of creation and 
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annihilation operators, 


= ^ / d^q Vat (g) Va,(g) (12) 

/ ^^9Vat(qr) Vaj^(g)+ Eo J d^q al{q) ay{q) (13) 

Winter = 9 j d^qal{q) {ay{q) + at(g)) a,,(qr) (14) 

with t denoting the adjoint operator, and ax^y{q) the annihilation operators for an x,j/- 
particle at location q in position space, formally defined by 



= Vm + l7p[{x^,q),y^) 

(15) 


= - g) \ Xi, y^) , 

(16) 

{ay{q)ip){x'^,y'^) 


(17) 

(4(9) 

1 

^ i=i 

(18) 


The combinatorial factors a/ttt. + 1 and arise from the fact that every conhguration 

of m x-particles occurs in m! different permutations. Since the operators a^^yiq) contain 
the point evaluation of functions, they can actually be defined rigorously on a dense 
subspace of while their adjoints are undefined for any nonzero V’- 

Model 1 is not Galilean covariant but a variant that is can be set up; it is defined, 
along with a version using IBCs, in ini- 


2.1.2 UV Divergence of the Original Hamiltonian 


The original Hamiltonian TTorig as defined in (10) above is UV divergent and thus ill 
dehned. The source of the difficulty is the delta function in the last line of (10): Since 


a delta function is not a square-integrable function, the right-hand side of (|10[) does not 
lie in for any choice of '0; thus, (10) does not dehne an operator in It can also be 


pointed out that the delta function 6^ (which plays the role in (10) of the wave function 
of a newly created y-particle) is a state of infinite energy. 


formally (^6^ 




2 m„ 




d-^k 


{hkf 

2m„ 


= oo 


(19) 


The standard procedure for obtaining a well-defined Hamiltonian (UV cut-off) is to 


7 


















replace the delta function by a square-integrable function <^9 : —>■ C, yielding 

fc 2 ^ fc 2 ^ 

^ i=l y j=l 

m „ 

+ c/\/n + 1 V / £yip*{Xi-y)7p{x"^,{y^,y)) 

1^1 

m n 

+ -%Y.Y1 - yj^ y" \ y^) ’ 

V i=l j = i 

which amounts to saying that an x-particle has an extended charge distribution with 
density function ip (if (p is real-valued). Equivalently, the interaction Hamiltonian is 
replaced by 

hfinter,cutoff — Q J d Q (l\.{(j) {(ly^ip{(j) ttxifl) j (^l) 

where 

(ay,UQ)^)(x"'yyl = +1 [ d^yip*{q-y)'ilj{x'^,{y'^,y)) (22) 

JR3 

1 

{al^{q)^){x"":y'') = -j=^^{yj-q)^{x^,y'"\yj ). (23) 

Vj = l 

2.2 Simplified Version: Model 2 

In Model 2, there is only one x-particle, and it is hxed at the origin. The only dynamical 
(and quantized) degrees of freedom reside in the ^/-particles; such models were considered 
by van Hove [HIE], and they tend to arise in the limit rux —>■ C) 0 . Specihcally, we take 
the conhguration space of Model 2 to be the set 

OO OO 

Q = Q„ = U 2»”’ = U (*" \ {“»” • P'l) 

n=0 n=0 

where 0 denotes the origin in a conhguration is denoted by y” = {yi ^..., y^). Note 
that we exclude the possibility that any ^/-particle can be at the location of the x-particle 
(i.e., the origin), but we do not exclude the possibility that two ^/-particles can be at 
the same location because for the purposes of this model, in which the ^/-particles do 
not interact, there is no need to exclude such conhgurations. 

Correspondingly, we take the Hilbert space to be the bosonic Fock space 

OO 

^ = 0 . (25) 

72=0 
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Since for the definition of the Hilbert space it plays no role whether the origin is excluded 
or not, elements xjj G can be regarded as comp lex-valued functions on Q that are 
permutation-symmetric in every sector; thus, C L^{Q), where Q is thought of as 
equipped with the measure /i dehned in the same way as in Q (L^(Q) contains also 
non-symmetric functions). The “original” Hamiltonian is now a simplihed version of 
0: 

fc2 

5 ^ Vyiv") + nEoiiy") 

y j=i 

+ gy/n+lilj{y^,0) 

n 

+ ^ ( 26 ) 

V 

Also this Hamiltonian is UV divergent. In the following, we will sometimes hnd it useful 
to write instead of '^(y"') in order to emphasize which sector of ip is being used. 


3 Interior—Boundary Condition and Corresponding 

Hamiltonian for the Two Models 

For Model 1 and Model 2, we describe the IBC and the Hamiltonian Hjbc- We also 
explain why Hjbc is a reasonable interpretation of the formula for i^orig- We begin with 
the simpler scenario of Model 2. 


3.1 IBC for Model 2 

Let denote the unit sphere in MP. The IBC demands the following: For every lo G 
ne{0,l,2 ,...}, 9“e(K3\{0})”, 


. (27) 

The IBC is a condition on the wave function pj at or near the “boundary” of conhguration 
space; the relevant boundary of is 


= 1J(M^)^“^ X {0} X . (28) 

i=i 


That is, the boundary consists of those conhgurations at which one of the ^/-particles 
collides with the x-particle. Due to the permutation symmetry of we can assume 

without loss of generality that it is the n + 1-st variable, Vn+i^ that approaches 0. The 
name “interior-boundary condition” reflects the fact that (27) relates the values (or 
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limits) of "0 on the boundary dQy to values of '0 in the interior of configuration space 
Q (namely, in the interior of a different sector corresponding to a lesser number of 
^/-particles). 

The IBC (27) allows in particular that ^ruj) diverges like 1/r as r —)• 0. In 


fact, it requires that V’ so diverges whenever the right-hand side of (27) is nonzero]^ By 
permutation symmetry, '?/’(?/"') diverges as any approaches 0. 

On wave functions 'ip satisfying the IBC (27), the Hamiltonian H = Hjbc is defined 
by 

fc2 ^ 


y j=i 


gVn +1 

4:71 


Im — ra;) j 






(29) 




This equation differs from the expression (26) for Tforig only in the second line, where 


0 ) has been replaced by a more complicated expression involving the behavior of 
Ip near the configuration (?/"', 0 ); after all, pj diverges at this configuration by virtue of 
the IBC, so the expression 'ip{y'^, 0) does not make sense. 


In view of the delta function appearing in (29), it may seem unlikely that such a 


Hamiltonian can be well defined. However, Theorem [T] in Section 3^ below shows that 
in fact it is well dehned. It turns out that the last line in (29) always gets canceled 


by contributions to the hrst line, which may contain delta functions because we now 
allow wave functions that diverge like 1/r = l/\yj\, and the Laplacian of 1/r is (as 
readers may remember from electrostatics, where 1 /r occurs as the Coulomb potential 
(p generated by a point charge, satisfying the Poisson equation A(p = dvrp with charge 
density p). As a consequence of this cancelation, Hjbc'P’ is a square-integrable function 
for pj satisfying the IBC. See also Remark in Section 3.3 below for further discussion. 

Let us point out a connection between the second line of (29), the line that dif¬ 
fered from the original Hamiltonian (26), and the corresponding line in (26). Think of 
' 0 (?/"',rn;) as a function of r; as r —)■ 0 , it can be expanded in the form pj{y'^‘ 
a/r + b + 0{r) with complex coefficients a, 6 ; then 


, rn; = 


d 

=b + 0 {r) , 


(30) 


^This divergence is to be expected if we keep in mind that j'bp represents probability density and note 
that for the inward radial flow in M®, i.e., for the dynamical system defined by the ODE x = —x/\x\^ 
the stationary density is 1/r^ = l/|a;p. To see this, note that if we squeeze a spherical shell of radius 
ri and thickness dr to a shell with smaller radius r 2 and equal thickness dr then its volume goes down 
(and thus, if its probability content is conserved, its density goes up) by a factor of [ri/r^)'^. That is 
why to expect \p\^ ^ 1/r^ and thus p ^ 1/r as r —>• 0. 
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which yields b in the limit r —0, and if '0 did not diverge as r —)■ 0, this would be 
exactly the value b = hmr\o t/’= t/’(2/”)0) occurring in (26) in the second line. 
So, (29) is really quite similar to (26). 

The Hamiltonian Hjbc is not the sum of two self-adjoint operators, the free Hamilo- 
nian and an interaction Hamiltonian. That is because the free Hamiltonian is defined 
on a different domain than Hibc^ containing wave functions that do not satisfy the IBC 
and do not diverge on the diagonal. It is the action of the Laplacian on functions that 
diverge on the diagonal that leads to delta functions and thus makes the last line of (29) 
possible, and it is only in conjunction with the IBC that the Hamiltonian (29) leads to 
conservation of probability. 

Indeed, in order to understand why the IBC (27) was chosen this way, and how it 
works together with the formula (29) for the Hamiltonian, it is illuminating to calculate 
the balance equation for and check that is conserved. We go through such a 
calculation in the next section. 


3.2 Self-Adjointness and Conservation of Probability 

Theorem 1. [9] On a certain dense subspace &ibc oj , the elements of which 

satisfy the IBC ( |^ , the operator Hjbc given by ( [2^ is well-defined and self-adjoint. 
If Eq > 0, then Hjbc positive. 


The relevance of this theorem is that it provides certainty that Hjbc is not afflicted 
by UV divergence. We give here a calculation checking on a non-rigorous level that 
probability is conserved. Using the symbol 


jy^ =3y^y'") = — 


(31) 


for the usual probability current, we obtain from (29) that, at any configuration y'^' 
(without any ?/-particle at the origin). 


d\ip{y^)\ 

dt 


X] ■ 3yj 


i=i 

gVn + l 
2Trh 


Im 




§2 


The last line can be re-written, using the IBC (27), as 


gy/nlH 27rh‘^\/n -\- 1 
Tm 


\im rfj* (y'^, ruj) I \im -^(r'l/jiy'^^ruj)) 

r\0 J r\o or \ / 


27ih 


g my r\o 

h 


d 


= — (n-|-1) lim / d^u) —ImrLj)—'0(|/"', ro?), 


r\0 




dr 


(32) 


(33) 

(34) 
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using that dr{r'ip) = 'ip + rdr'i/j and \m.{r\'ip\^) = 0. Thus, 

d\Hy'')l_ _ _ ^ ^ 1 ) ^2 f ,2 

y3 ■•Vj ^ ^^^0 

§2 


dt 


y, ■ 3y^ - (n + 1) lim r J d u; u; ■ {y^, ru)). (35) 


This equation possesses a simple interpretation: —times the integral is just the flux 
of probability current toward the origin across the sphere of radius r in the coordinate 
space of Vn+ii f^e limit of that as r —)■ 0 is the current of probability into the origin 
in the coordinate space of Hn+i] summing over all ^/-particles would yield, due to the 
bosonic symmetry of "0, n + 1 equal terms; thus, the second summand on the right-hand 


side of (35) is the total flux of probability into the boundary configurations obtained 
from r/” by adding one ^/-particle at the origin. That is, the balance equation (35) asserts 


that the probability density |'0r changes in two ways, due to transport of probability in 
the n-particle sector and by increasing at just the rate at which probability disappears 
on the {n + l)-particle sector by flowing into the boundary dQy^\ Therefore, total 
probability is conserved. 


3.3 Remarks 


1. Comparison to Bethe-Peierls boundary condition. The IBC (46) has some parallels 


to the Bethe-Peierls boundary condition [2], which reads for a wave function ih e 

d 

^ ^rip{r(jj'^ = 0 (36) 


r ( 9 

hm a + ^ 
r\o V or 


with given constant a G M. This condition is used for giving precise meaning to a 
Schrodinger equation for -0 : —)■ C with a Dirac delta function as the potential. 


H = -£;& + gS\x), 


(37) 


a kind of interaction known as a point interaction [T]. Like our IBCs (27), 


(46), the Bethe-Peierls boundary condition (36) concerns the “boundary” at r = 0 


However, in contrast to the IBC, which connects two sectors of Q, the Bethe- 
Peierls boundary condition involves only one sector, as the wave function -0 is 
defined on Also, the Bethe-Peierls boundary condition implies zero current 
into r = 0 (see, e.g., [16] for more detail), whereas the IBC leads to a nonzero 
current into the boundary. 


2 . Comparison to vertex conditions in networks. An example of a known boundary 
condition leading to nonzero current into the boundary is provided by the ver¬ 
tex conditions for quantum mechanics on networks (see, e.g., [H] and references 
therein). A network, or graph, consists of several 1-dimensional spaces (i.e., in¬ 
tervals) called edges, glued together at their end points called vertices. The wave 
function is a (say, complex-valued) function on the network (i.e., on the union of 
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the edges), and the Hamiltonian is given by the Laplace operator on each edge, 
along with bonndary conditions for the end points of the edges, also called vertex 
conditions. The vertex conditions are related to the fact that, since no positive 
amonnt of probability can be located at a vertex, all of the probability cnrrent 
into the vertex must be compensated by a current out of the vertex; that is, for 
each vertex n, the sum of the currents along all edges e connected to v (pointing 
away from v) must vanish, 

(38) 

e 

While Dirichlet or Neumann boundary condition at every end point of every edge 
would lead to jeiy) = 0 for all v and e, it is also possible to impose vertex conditions 
on 'ifj that allow nonzero flow of probability from one edge across a vertex to another 
edge, the simplest one being the conjunction of 


lim 'ip{x) = lim 

along e x^v along e' 


(39) 


i.e., that -0 is continuous at vertices, and 


d,iiv) = aiiv), 


(40) 


where a G M is a given constant and de'ip{v) means the derivative of V’ along the 
edge e, taken at v. 


To draw parallels between (40) and IBCs, we may compare the edges of a network 


to the sectors of the conhguration spaces considered in this paper; let us call both 
the “blocks” of the space. Of course, the edges of a graph have equal dimension 
while the sectors of our Q have different dimension. We may note that in both 
situations the probability loss in one block is compensated by a probability gain in 


others. Also, the condition of continuity (39) is a relation between the boundary 


value in one block and a value in another, but in contrast to an IBC the latter is 
also a boundary value; that is, in analogy to the terminology “interior-boundary 
condition,” (39) is a “boundary-boundary condition.” 


Combinatorial factors. As remarked at the end of Section 2 .1.1 the combinatorial 
factors ^/n and \/n + 1 that appear in the Hamiltonian arise from the fact that 
we use ordered conhgurations ..., although in nature conhgurations are 
unordered, as in {y ^,..., y^}. This is also true of the IBC (27) and the associated 


Hamiltonian (29). If we used unordered conhgurations, which leads to topologi¬ 
cally non-trivial conhguration spaces (see, e.g., mm, these factors would not 
appear. 


4. Neumann-type IBC. Two well known types of boundary conditions are Dirichlet 
boundary conditions. 




dQ 


= 0 , 


(41) 
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and Neumann boundary conditions, 


dtp 

dn 


dQ 


= 0 , 


(42) 


with n the normal vector on the boundary hypersurface. While the IBC (27) 


resembles the Dirichlet type in that it involves the values of ip (or rather, rip) on 
the boundary, one can as well set up a different IBC that resembles the Neumann 
type in that it involves the normal derivative of ip (or rather, rip). Here, r = 0 
plays the role of the boundary, and the radial direction plays the role of the normal 
to the boundary. The Neumann-type IBC asserts that for any n G {0,1,2,...}, 
any conhguration ?/"' G (M^ \ {0})"', and any uj G 




ruj] = 


gmy 




2Trh?\/n + 1 

Also this IBC typically leads to ip{y"',ru:) diverging like 1/r as r —)• 0. 


(43) 


On wave functions ip satisfying the IBC (43), the Hamiltonian H = Hjbc is defined 
by 


{Hip)iyn 


2m 


■,p(y’') + nEa;P(y") 


y 7=1 



(44) 


This equation differs from (29) (for the Dirichlet case) in the second line, which 


no longer involves a derivative, and in the last line, which we will discuss further 
in Remark [6] below. 


While (43) and (44) together define a different time evolution for ip than (27) and 


(29), Equation (35) is still true and guarantees that the amount of probability lost 


on the n + 1-sector of Q due to probability flux into the boundary is added on the 


n-sector. This fact suggests that also (43) and (44) dehne a self-adjoint operator 


5. Robin-type IBC. A Robin boundary condition is one of the form 




(45) 


with given constants a and (I. An IBC of an analogous form can be set up as 
follows: for any n G {0,1,2,...}, any configuration y'^ = {y-y, ..., y^) G (M^\{0})” 
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of ^/-particles, and any uj G 


lim ( a + (3^ ) ( rijj{y'^, ru;) ) = 


d 


r\0 y 


dr 


2m. 


h^y/n + 1 


i’iy") . 


(46) 


where {a, (3) G \ {(0,0)} are constants. Also this condition typically leads to 


^jJ that diverge like 1/r as r —)■ 0. The Dirichlet-type condition (27) is included 


in this scheme for a = —A'K/g,[3 = 0, while the Neumann-type condition (43) is 
included for a = 0 , /9 = dvr/gf. 

The associated Hamiltonian H = Hjbc is given by 

fc 2 ^ 


2m 


y j=i 


(fuj lim^7-|-(5—j rcj) j 

S2 

n „ 

\ ruj) , (47) 

.7 = 1 L 


2m. 


where the constants 7 , 5 G C satisfy 


a*7 G I 
P*5 G I 
a*5 - 7*/3 = 


= - 1 . 


(48) 

(49) 

(50) 


It should always be obvious when the symbol 5 means the constant 5 G C and 


when the Dirac delta function. Note that for 7 = 0, 5 = the expression (47) 


for the Hamiltonian agrees with (29), while for 7 = gj^r:^ <5 = 0 it reduces to (44) 


On a non-rigorous level, the conservation of probability (i.e., self-adjointness of the 


Hamiltonian (47)) can be checked by means of a calculation similar to the one in 
Section |3^ above and to the one in Section 4^ below. The self-adjointness breaks 
down (already on the non-rigorous level) if we relax the conditions (48)-(50). 


6 . Dirac delta function terms in H. Here is a reason for thinking that, among the 
many different IBCs that are mathematically possible corresponding to different 


choices of the constants a,(3,'j,S, only the Dirichlet-type IBC (27), corresponding 


to a = —An/g, /? = 0 = 7 , and 5 = g/An, is physically relevant as a replacement 


of the original (UV divergent) Hamiltonian (26): It is the only choice that leads 


to a term in Hjbc that reproduces the Dirac delta function terms in Hor\(r, i.e 


for which the last line of (47) agrees with the last line of (26) 


orig5 

That is because 

the Dirichlet case is the only case in which the last line of (47) can be expressed 
in terms of _ 
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And the last line of (47) is dictated by the condition that Hip has to be an H 


fnnction and thns cannot contain contribntions that are Dirac delta fnnctions. 
Indeed, since for any a,(3,'y,5, the wave fnnction ip^'^'> diverges at the bonndary 
like 1/r, the Laplacian in the Hamiltonian always yields a distribntion of the form 

- E \ + siy") (81) 

j j 

with some fnnctions / and g] so the Dirac delta contributions need to be canceled, 


which leads to the last line of (47). 


7. Positivity. Another reason for thinking that the Dirichlet-type IBC (rather than, 
say, Neumann-type) is the physically relevant choice may be that Hjbc for a 
Dirichlet-type IBC is positive, as mentioned in Theoremj^in Section 3^ Generally 
speaking, the IBC approach neither requires nor guarantees that Hamiltonians 
are bounded from below. Presently, we do not know for which other choices of 
a, f3, 7 ,6 the Hamiltonian will be positive, but we see reason to believe that, in the 
variant of the equations appropriate for Model 1, the Neumann-type IBC leads 
to a Hamiltonian that is not bounded from below. Further considerations about 


physical reasonableness of Hjbc can be found in Sections 3.5 and 3.6 below. 


8 . Bohmian trajectories. There is a natural way of dehning Bohmian trajectories for 
the models described in this paper; we describe this in detail elsewhere [5]. The 
Bohmian configuration Qt follows a Markov jump process in configuration space 
that is IV’tp-distributed at every time t. The process has finitely many jumps 
in every finite time interval and moves deterministically in-between according to 
Bohm’s equation of motion. In our Models 1 and 2, the jumps correspond to 
the creation or annihilation of a particle. The jumps to a lower sector (particle 
annihilation) occur whenever a ^/-particle hits an x-particle; in that event, the y- 
particle gets deleted from Qj. While the jumps to a lower sector are deterministic, 
the jumps to a higher sector (particle creation) are stochastic. They can occur 
at any conhguration Qt, with a rate depending on Qt and ipt, and lead to a 
conhguration with a new y-particle created at the location of an x-particle; the 
^/-particle then moves in a random direction that is uniformly distributed over 
the sphere. Compared to previous models of particle creation and annihilation in 
Bohmian mechanics (see [18] and references therein) that involved a UV cut-off, 
the difference is that in Models 1 and 2, the ^/-particle gets created at (rather than 
near) an a:-particle, and that annihilation in deterministic. The process is time- 
reversal invariant (notwithstanding that annihilation is deterministic and creation 
is stochastic). 


3.4 IBC for Model 1 

We now describe an IBC and the corresponding Hamiltonian for Model 1. The IBC 
demands that for any m,n E {1, 2 ,...}, any configuration = {xj ,..., Xm) E of 
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x-particles, any configuration , y^) G of y-particles with fl j/"' = 0 

(i.e., Xi 7 ^ y^ for all i, j), any i = 1 ,..., m, and any j = 1,..., n, 

lim \y^ - Xi\'ilj{x'^, y^) = On-i ^J^Xi = x, y ^), (52) 


where ^ denotes omission, 


Q^n —1 


m^ruy 


(63) 


2'Kh?y/nmx + my ’ 

and G M is the same coupling constant as before. 

The IBC is a condition on the wave function -0 near the diagonal A in configuration 
space Qx x i.e., the set of “collision conhgurations,” 


A = 1(1’", y") € Q, X Qy : Xi = y, for some i. j J . 


(54) 


If we regard the collision conhgurations in A as not admissible conhgurations then the 
conhguration space is the set diherence Q = (Qx x Qy) \ A, and its “boundary” is 
dQ = A. 

On wave functions -0 satisfying the IBC (52), the Hamiltonian H = Hjbc is dehned 


by 


(HiBci^)(x^.y^) = 

+ nEo'ilj(x^,y'^) 

, gVnTT 




2m, 


E 9") - ^ 9") 


2=1 


y j=i 


dvr 


/ d^u) lim — r'tplxi Xt - y^^^ -^Xi + yyVu'^ 

^=ll, r 0 ry 




i=l j=l 


with 


m„ 


yx 


mx + m. 


and jjLy = 


m^ 


mx + m. 


(55) 


( 66 ) 


We conjecture (and are working on a proof [TO]) that Hjbc given by (55) dehnes a 
self-adjoint operator on a suitable dense subspace of ® , the elements of 

which satisfy the IBC (52). 


3.5 Ground State Energy and Effective Yukawa Potential 


Elsewhere [9| we show that Hjbc as in (29) for Model 2 with Dirichlet-type IBC (27) 
possesses a non-degenerate ground state V'min- It is given by 

<*..,„(9......y„)=JV^(-^) n -rn- (V) 


i=i 


12/7 
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with normalization constant 


AT = exp 


47rh^ ^j2myEQ / 


(58) 


and eigenvalne 


g^my^J2myEQ 

2'Kh^ 


(59) 


That is, the state is a snperposition of different nnmbers of ^/-particles, and in each sector 
all ^/-particles have the same wave fnnction; so the x-particle at the origin is dressed 
with a cloud of ^/-particles. The probability distribution of the number n of ^/-particles 
is a Poisson distribution with mean value 


(n)p. 


g^ml 

2TThA^2myEQ 


(60) 


We claim further that in Model 1, in which the x-particles interact by exchanging 
^/-particles, the a:-particles effectively interact through a Yukawa potential; this in fact 
agrees with the result of Yukawa’s original reasoning [20], see below. A simple way of 
computing the effective interaction potential is to consider (as “Model 2b”) N x-particles 
hxed at a?!,..., G and to hud the ground state energy. To this end, consider wave 
functions of j/-con£gurations that simultaneously satisfy N IBCs, 


lim ( 


r\0 \ 


Xi -h rn;) = 




2TTh^^/n + 1 




(61) 


for every i = 1,... ,N, u> G n G {0,1, 2,...}, and y” G (M^ \ {a^i,..., a^w})"'- The 
corresponding Hamiltonian reads 

fc 2 ^ 


y j=i 


gVn + l 


Att 


N . ^ 

/ d^u) \im.— ^Xi + ru 

“ J r\o dr \ 

S 2 


N n 




( 62 ) 


i=l j=l 


For Eq > 0, the ground state is 


n N ^-^2myEo\yj-Xi\/h 

'IprainiVl,- ■ ■ ,yn) = CnYlYl -^-Z-^ 

j=i i=i \yj 


18 















with suitable factors and eigenvalue 


Err\\r\ — 


g^rriy 


/ N^j2myEo 


- E 


^-^2myEo\xi-Xj\/h 


\Xj — Xn 


(64) 


\ 2h 

That is, the ground state energy of the ^/-particles, given the ^(-particles at Xi,... ,xn, 


is given by (64). Regarding this energy function of a^i,..., a^Tv as an effective potential 


for the x-particles (which is appropriate when the ^(-particles move slowly), we see that 
^(-particles effectively interact through an attractive Yukawa pair potential, 

g—A /2 

V{R) = const.-— (65) 

R 

with R the distance between two x-particles. If we take the energy needed to create a 
^/-particle to be Eq = niyC^, then 


A = 72^ 
n 


( 66 ) 


which is, up to the factor a/2, the value originally obtained by Yukawa [20] considering 
the effective interaction of nucleons by exchange of pions. We expect that the factor a/ 2 
is owed to the non-relativistic nature of our model. 


3.6 IBC Hamiltonians as a Limit of Removing the Cut-Off 


If one introduces a UV cut-off into the UV divergent original Hamiltonian (26) of Model 
2 , it becomes the well-dehned operator Hy, given by 


HMvl = 


y j=i 


2 m, 

-h gVn + 1 


v’{y)i’{y’',y) 


9 




(67) 


i=i 


analogous to (20) for Model 1. Here, the Dirac delta function 5^ has been replaced by 


(p : —)■ C, a square-integrable function describing the charge density of the x-particle, 
and the limit ip ^ would correspond to removing the UV cut-off. It is well known 
uniEi that, if Eq ^ 0, there are numbers Ey^ ^ 1^ such that the operator Hy^ — Ey^ 
possesses a limit as —)■ 6^; let us call this limit, which is a self-adjoint operator. Hoc- 
Since Ey, —)■ cx) in this limit, Hy, tends to inhnity in a sense, but this sense is harmless 
because two Hamiltonians that differ only by a multiple of the identity operator can 
be regarded as equivalent, as they generate the same time evolution (if we regard wave 
functions that differ only by a global phase factor as equivalent). So, in a relevant 
sense, the limit ip ^ can indeed be taken, and Hoo is the limiting Hamiltonian, which 
suggests regarding Hoo as the “physically correct” Hamiltonian. 
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Theorem 2. [H] For the Hamiltonian Hjbc of the Dirichlet type for Model 2, defined 

Hoo = Hibc + E^o ( 68 ) 


hy (27) and 


with the constant E^o = ^/Eq/4:7^ ■ 


On the one hand, this result lends further support to regarding Hibc and the IBC 
approach as physically reasonable. On the other hand, the result provides a more direct 
representation of H^o than was available so far. 


4 IBCs For Codimension-1 Boundaries 

In Model 1 and Model 2, our examples for the application of IBCs to QFT, the boundary 
of each sector of conhguration space was a subset of codimension 3, corresponding to 
an x-particle and a ^/-particle being at the same location. IBCs can also naturally be 
used for codimension-1 boundaries in conhguration space, as we outline in this section. 
A simple example we describe elsewhere [16]. As far as we can see, the considerations in 
this section have no impact on avoiding the UV divergence without UV cut-off, but they 
form a natural framework, although not the most general one, for IBCs. Preliminary 
considerations in this direction were described in [TH]. 

4.1 Configuration Space and Hilbert Space 

We take the conhguration space Q to be a hnite or countable union of disjoint manifolds 
with boundary, Q = (By dehnition, in a manifold with boundary, a neighbor¬ 

hood of an interior point looks like a piece of while a neighborhood of a boundary 
point looks like a piece of a half-space in In particular, the boundary has codimen¬ 
sion 1, i.e., dimension d — 1. The boundary may be empty.) We write for the 

boundary of dQ = UndQ^'^f and Q° = Q \ dQ for the interior of Q. We take Q 
to be equipped with a Riemann metric gij, which also dehnes a volume measure on 
and thus a measure jj, on Q, /i(S') = likewise, the metric dehnes 

a surface area measure \ on dQ. 

Wave functions can be complex-valued functions on Q. However, we can also be 
more general and allow cross-sections of vector bundles. Readers unfamiliar with vector 
bundles may ignore this further generality and think of complex-valued wave functions. 
So, for every n, let be a Hermitian vector bundle over of hnite rank = 

dime Eq"''^ (dimension of the hber spaces), i.e., a complex vector bundle equipped with 
a positive dehnite Hermitian inner product ( , )ij in every hber Eg^\ q G and a 

metric connection (i.e., a connection relative to which the inner product is parallel, or, 
equivalently, a connection such that the parallel transport it dehnes along any path from 
q to q' is a unitary isomorphism E^'^ —)■ We write E for and Eg for 

if q E The wave function will be a cross-section of E, i.e., a mapping %Ij ■. Q ^ E 

such that fj^q) G Eg for every q E Q. 
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The Hilbert space M’ = L'^{Q, E, consists of the square-integrable cross-sections 
of E and is eqnipped with the inner prodnct 

(^10) = [ Kdq) (69) 

JQ 

Note that Jq means the same as ~ 


4.2 IBC and Hamiltonian 


The IBC will be so constructed that the amount of probability per time that flows out 
of the boundary at q' E dQ gets added to IV'P at an interior point 

q = f{q') (70) 


in a different sector, / : dQ —j- Q°. We suppose that the derivative of / has full rank, 
i.e., that the image of / in does not have lower dimension than in particular, 
if q' E and q = f{q') E then dim < dimcIQ‘^”'0 = dim _ x_ Since 

many boundary points q' can be mapped to the same interior point q, the set of which 
will be denoted 

f~^{q) = {q' e dQ ■. f{q) = q }, (7i) 

we will need to make use of a measure over f~^{q)- The appropriate (unnormalized) 
measure for our purpose is 


z/q(-) = weak-lim 


A(- nf-\dq)) 
dq^{q} P-{dq) 


(72) 


For example, if S' = / ^(g) is a submanifold of dQ^'^'^ of dimension k and dQ^'^'^ 

has dimension £, then the density of Uq relative to the volume measure arising from the 
Riemann metric on S is 


det (yQQQ(n') (Cj, Cj)^ 


,3<t- 


det ( 5 ( 5 ( 6 *, Cj 




det (gQ(n) {df{ei),df{ej)) 




1/2 


(73) 


where df : Tq/dQ -E TqQ is the derivative of /, and e* are any linearly independent 
vectors in T^dQ with the hrst k in Tq/S] the quantity (73) does not depend on the 
choice of 6 *]^ In particular, if f~^{q) is a hnite or countable set (say, / is a local 
diffeomorphism), then for any q' G f~^{q), 

- 1/2 


^q{{q'}) = ^ lim 

rin '—J n' 


\{dq') 


dq'^{q'} nifidq')) 


det[gQ(r.) {df{ei),df{ej)) 


i,j<^ 


(74) 


^Alternatively, i^q can be expressed as a differential form Vq of maximal degree on S, 
Vq{q'){vi,... ,Vk) = \{q'){vi,... ,Vk,ei,... ,et) for any q' E S and vi,...,Vk G Tq>S, where A is the 
differential form corresponding to the measure A (i.e., the Riemannian volume form on i 9QF and 
the Ci are any vectors such that df(ei ),..., df{ei) is an orthonormal basis of TqQ. 
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for any orthonormal basis {ej} of the tangent space Tg/dQ. 

We now set np the Hamiltonian and IBC. We may inclnde a potential, either as a 
fnnction 1/ : Q —)■ M or more generally as a cross-section oi E ® E* that is pointwise 
self-adjoint; here, E* denotes the dnal space of Eg, so an element V{q) of Eg ® E* 
corresponds to an endomorphism Eg —)■ Eg that we also denote by V{q), and V being 
pointwise self-adjoint means that V{q) is self-adjoint on Eg relative to ( , )g. 

The IBC demands that for every bonndary point q, 

(a{q) + /3{q)dn'^i^{q) = p^(/(g)) , (75) 

where a{q) and (3{q) are given complex-linear mappings Eg —-E/(q) and means the 
normal derivative, i.e., the directional covariant derivative in the inward normal direction 
to the bonndary at q (normal in terms of the Riemann metric Qij). 

The Hamiltonian is, for any interior point q: 

Hil){q) = -^Aip{q) + V(q)^(q) + f i^qidq') (^j{q') + S{q')^r^^^p{q') . (76) 

f-Hg) 

Here, A is the Laplace operator associated with the Riemannian metric of Q and the 
connection of E (see, e.g., [1] for a detailed dehnition), and the coefficients 7 (g') and 
6{q') are given complex-linear mappings Egi —)■ (The coefficient 6{q') should not 

be confused with a Dirac delta function.) The functions a, /3, 7 , 6 are required to satisfy 


at every q G dQ the conditions 

Q;(g)^ 7 (g) : Eg ^ Eg is self-adjoint relative to ( , )g (77) 

/?(g)^5(g) : Eg ^ Eg is self-adjoint relative to ( , )q (78) 

a{qy 6{q)--f{qy (3{q) =-Ie,, (79) 

where Ie^ means the identity operator on Eg and a^qY : 77/(q) —t Eg means the adjoint 
of a{q) relative to ( , )q and ( , )/(g), i.e., 

(x, a(g)V)g = (a(g)x, 4>) j(^g) (80) 

for any x ^ Eq and 0 G 77/(g); likewise for /3(g)7 

We think of the masses as incorporated into the metric Qij, as in, e.g., 

= rriidxl + rriidyl -|- rriidzl + m2dx\ + m2dy2 + m2dz2 (81) 


for two particles of different mass in Euclidean space (see jl| for further discussion). 
Then the mass need not be put into the prefactor of the Laplacian (in the Hamiltonian 
77 as in (76) above) or the gradient (in the current j as in dSd]) below) 


^This convention has the possibly undesirable consequence that, when different sectors correspond 
to different particle number, the Riemannian volume y is weighted in different sectors with different 
powers of the mass (such as m^"); however, this can easily be compensated by reweighting by a 
factor of which in turn requires, if f{q) contains one particle less than q, a further factor of 

J71-3/2 a,l3,^,S. 
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This completes the dehnition of the Hamiltonian. In contrast to Model 1 and Model 
2 with IBCs, 'ijj does not diverge at the boundary in this case. 


4.3 Conservation of Probability 

Here is a formal (non-rigorous) derivation of the conservation of i.e., a check of self¬ 
adjointness of H on the non-rigorous level. By we mean 

which is the density relative to fi of the probability distribution in Q associated with 
Ip G with 11-011 = 1. It evolves in general according to the balance equation 


= 2 lm{^p{q), H^/J{q))^ . (82) 

It is known (e.g., |1]) that for H = ^ Hermitian vector bundle over a Riemannian 

manifold, 

I lm{^p{q), H^/J{q))^ = -div j{q) (83) 

with the probability current vector field 

j{q) = him (V^(g), VV'(g))^ (84) 


on Q. Here, divj denotes the divergence of the vector held j; in coordinates, divj = 
Xla where Da is the covariant derivative operator arising from the Riemann metric 
on Q. The gradient Vp! is the T^-valued vector held obtained from the T^-valued 1-form 
that is the covariant derivative of pj by “raising the index” using the Riemann metric. 


Now for the Hamiltonian (76), the balance equation (82) becomes 


dt 


= -divj(g) + / ^g{dq') I Im(V'(g), [7(g') + 5{q')dn]'ip{q') 


(85) 


f-dq) 


For each g', the integrand can be written, by virtue of the IBC (75), as 


l\m{^^[a{q) + I3{cl)dn]p){q), [l{cl) + d{q)dn]p){q))^ 

= hIm(^?/;(g'),a(g')W(?')^(?'))^, 

+ h\m{^{q'),a{q')^{q')dn'ip{q')^ ^ 

- hlm(^^p{q),-i{qyP{q')dn'ip{q')^ 


( 86 ) 


(87) 


where the hrst and the last term vanish by virtue of (|73 and ([78|). What remains is 


hlmU{q), a{q)^5{q)--i{q)^(5{q) dn'ipiq) 


q' 


( 88 ) 
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which agrees, by virtue of (79), with 


— him 


/), ))^, = -jn{q '), (89) 

where jn{q') means the component of j{q') normal to the boundary, or 

jn = f Qij (90) 

with n the inward-pointing unit normal vector to the boundary. Thus, in total, 

d\i){q)\^ 


dt 


(91) 


=-divj(g) - / ^q{dq) ]n{cl) ■ 

Now, if in{q') < 0 then —jniq') X{dq') dt is the amount of weight lost in the sector 
containing q' due to current into the boundary region dq' around q' within duration 
dt. Likewise, if jn{q') > 0 then jn{q') X{dq') dt is the amount of weight gained in 
the sector containing q' due to current coming from dq' within duration dt. That is, 
jn{q') X{dq') dt is the net gain, positive or negative. Now the second term on the right- 


hand side of (91) represents a gain in the amount of weight (while the divj term 
represents transport of weight within one sector); in fact, the gain in the region dq 
around q within duration dt is 


-li{dq)dt / Uqidq) jn{q') =-dt / X{dq') jn{q). 


(92) 


f-Hg) f-Hdq) 

Thus, the net gain in dq exactly compensates the net loss in f~^{dq), and 
Jq lV’('?)Ph('^'?) is conserved. 

Equation (91) can be regarded as a transport equation for the weight, with 
two types of transport: continuous motion within a sector of Q, and transport between 


sectors of Q (either from q' to f{q') or from f{q') to q'). Equation (91) actually is a 
probability transport equation for the |'0p-distributed stochastic process in Q described 


m 


It also becomes clear from the above derivation of (91) that the conditions (77), 


(78), and (79) cannot be weakened within our scheme without losing (91) and thus 


the self-adjointness of H. (Except that (77)-(79) may fail on a A-null set of boundary 
configurations, or a,/3,'y,S may be undefined on such a set.) After all, if gpi^q') and 


dn'>p{q') can be chosen independently, then the only way in which ( 86 ) can always be 


equal to (89) is if (77), (78), and (79) are true. Now V’(<?0 ^^cl dni^{q') can be chosen 
independently by appropriate choice of initial data for —despite the IBC (75), which 
can be satisfied by appropriate choice of 'ip{f{q')). To be sure, (91) can be true for all 


"0 satisfying the IBC also if the integrals in (85) and (91) agree while the integrands 
are not equal. For example, this happens when 7 = 0 = 5 (so (79) is violated), f~^{q) 
contains two boundary points, say q' and q", and the loss at q' always compensates the 
gain at q” and vice versa (e.g., if jniq') = —jniq") and i^qi{q'}) = ^qi{q"}))- However, 
such possibilities lie outside our scheme, according to which the weight lost at q' is added 
to fiq'), and will not be considered here. 
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4.4 Application: Cut-Off Radius 

IBCs on a codimension-1 boundary can be used for implementing an unusual kind of 
UV cut-off, in which the x-particle is smeared out, but not over a ball but over a sphere 
of small radius 5 > 0. While an a:-particle smeared out over a ball can emit and absorb 
^/-particles anywhere within that ball, an x-particle smeared out over a sphere can only 
emit and absorb ^/-particles at a distance from its center that is exactly 6. That is, a 
^/-particle gets absorbed as soon as it reaches distance 6 from the center of an x-particle, 
and gets created on the sphere of radius 6. We exclude the possibility that any ^/-particle 
is ever closer than 6 to an x-particle. This kind of UV cut-off was hrst described, as far 
as we know, in [TH]; we will call it a “5-cut-off” in the following. 

When we apply this cut-off to Model 1, the conhguration space is 

Q = e Q^x Qy-. \xi - yj\ > 5 Vi, j| (93) 

with Qj. and Qy as in ([^, and its boundary 

dQ = e Q : \Xi — yj\ = 5 for some i, j| (94) 

has codimension 1 almost everywhere. For Model 2, the conhguration space would be 

Q={l/"eQ,:|?/,|>5Vj} (95) 

with boundary 

dQ = I?/”' G Qy : \yj \ = 6 for some j| (96) 

(also of codimension 1), and the Hilbert space is = L‘^{Q). Let Bs denote the 5-ball 
around the origin, = {^ G : \y\ < 5}. 

The IBC (for simplicity for Model 2 with 5-cut-off) of Dirichlet type demands the 
following: For every cj G n G {0,1, 2,...}, y"' G (M^ \ Bs)^, 

iTTh^oQn + 1 

The associated Hamiltonian is dehned by 




2m 


y j=i 


gy/n + 1 
dvr 


,2 d 
or 


-=5 




( 98 ) 


at any y"' e Q \ dQ. 
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In the language of Sections 4.1-4.3 


9ij 

vector bundle E is the trivial rank-1 bundle E 


iTbydij, f{y"‘) = y"' \ Bs, the Hermitian 
= Q X C, yU is the volume as in ([^, 


= 1J(M3 \ X §2 X (M3 y ^ 

i=i 


A is locally vol3(j_i) x area x vol3(„_j), 


/ Su;, Vj, • • •,?/n); ^ e , (lOO) 

i=i 


Vyri^dPuj) = 52 (ptjj on any of the n spheres in (100), «(?/”') 
0, 7(2/"') = 0, and 5{y'^) = gmy/A7i5\/n + 1. 


-dvr^A/n l/gmy, (5{y'^) = 
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